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Based on the color-spin interaction in diquarks, we argue that charmed multiquark hadrons are 
likely to exist. Because of the appreciable number of charm quarks produced in centrál nucleus- 
nucleus collisions at ultrarelativistic energies, production of charmed multiquark hadrons is expected 
to be enhanced in these collisions. Using both the quark coalescence model and the statistical 
hadronization model, we estimate the yield of charmed tetraquark meson Tec and pentaquark baryon 
Ocs in heavy ion collisions at RHIC and LHC. We further discuss the decay modes of these charmed 
exotic hadrons in order to facilitate their detections in experiments. 
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I. INTRODUCTION 

Possible existence of exotic mesons consisting of two 
quarks and two anti-quarks was first suggested by JafFe 
in the framework of the MIT bag model [I . Since 
then, there have been continuous discussions on whether 
the mesons in the scalar nonet are candidates for such 
tetraquark mesons. Recently, interest in tetraquark 
mesons has been extended to include those containing 
heavy quarks as several heavy mesons, that were 

observed in B meson decays, do not seem to fit well 
within the conventional quark model |3]. Tetraquark 
mesons with two heavy anti-quarks {QQqq), henceforth 
called Tqq , are particularly interesting as they are explic- 
itly exotic from flavor considerations [5]. Moreover, sim- 
ple theoretical consideration based on the color-spin in- 
teraction [6J shows that for such configurations the bind- 
ing energy increases as the mass of the heavy quark in- 
creases. Calculations based on the flavor-spin interaction 
[7J Hl IS] O'" the instanton induced interactions [TD] also 
show that the mass of Tec is below that of two charmed 
mesons. For a similar reason, the chance of having a sta- 
ble heavy pentaquark (qqqqQ) increases as the mass of 
heavy anti-quark becomes larger. 

Experimental observation of such explicitly exotic 
hadrons is crucial in refining our understanding of mul- 
tiquark interactions in low energy QCD. However, pro- 
ducing the Tqq from an elementary process is highly 
suppressed as it involves creating two QQ pairs from 
the vacuum. In contrast, in relativistic heavy ion col- 
lisions at LHC, čc pairs are expected to be abundantly 
produced Since the hadronization from the quark- 
gluon plasma produced in these collisions tends to foUow 
a statistical description, production of exotic hadrons in 
heavy ion collisions at LHC is thus much more favorable 
than in elementary reactions [12l [131 E] • 
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In this work, we first give a qualitative argument why 
multiquark hadrons consisting of heavy quarks are likely 
to exist. Using both the quark coalescence model and the 
statistical hadronization model, we then give estimates of 
how many Tqq and charmed pentaquark baryons, if they 
exist, will be produced in centrál heavy ion collisions at 
both RHIC and LHC. Furthermore, possible decay modes 
of these charmed exotic hadrons are discussed. 



II. A SCHEMATIC MODEL FOR HADRON 
MASS DIFFERENCES 

A. known hadrons 

Sophisticated constituent quark model calculations 
have been performed to study possible stable multiquark 
hadrons that consist of heavy quarks. These results can 
be roughly understood in terms of simple arguments 
based on the color-spin interaction. To illustrate the 
mechanism, we introduce the foUowing simplified form 
for the color-spin interaction [6]: 

CHY.S^■S-—. (1) 

l>j J 

Here m and s are the mass and spin of the constituent 
quarks i and j. The strength of the color-spin interac- 
tion Ch should depend on the wave function and the 
exact form of the interaction as well as the color struc- 
ture of either the quark-quark or quark-antiquark pair. 
The color factor would be 8/3 for diquarks in the color 
antitriplet channel and 16/3 for quark and anti-quark 
pair in the color singlet channel. This simple form with 
Ch = Cb for a diquark and Ch = Cm for a quark- 
antiquark pair can capture some of the essential physics 
in hadron masses. To illustrate this point, we assume the 
foUowing constituent quark masses: niu^d — 300 MeV, 
rUs = 500 MeV, = 1500 MeV, and m 'b = 4700 MeV. 

Table |l] shows the mass differences between baryons 
that are sensitive to the color-spin interaction only. By 
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Mass DifT. Ma - Mjv My. - Ma Me, - Ma^ M^^ - Ma^ 
Formula |% % (i _ Hl- ) %(i _ 221" ) % (i _ zhil) 

Fit 290 MeV 77 MeV 154 MeV 180 MeV 
Experiment 290 MeV 75 MeV 170 MeV 192 MeV 

TABLE I: Baryon mass relations. The first column is fit to 
experiments. 

fitting Cb to Ma - M^y, we obtain C^/m^ = 193 
MeV and find that the mass difFerences Ms — Ma and 
Ms^ — Ma^ are well reproduced. This is in no way an 
attempt to make a best fit, but the point is that with 
typically acccpted constituent quark masses, the mass 
sphtting larger than Cb, refiecting that the quark and 
anti-quark correlation is about 3 times stronger than that 
between two quarks. 



Mass Diff. 


Mp - IvU 


Mk' - Mk 


Md' - Md 


Mb* - Mb 


Formula 
Fit 

Experiment 


Cm 

635 MeV 
635 MeV 


Cm 




Cm 


381 MeV 
397 MeV 


127 MeV 
137 MeV 


41 MeV 
46 MeV 



TABLE IL Meson mass relations. The first column is fit to 
experiments. 



When both quarks are heavy, the value of Ch is ex- 
pected to become larger as the strength of the relative 
wave function at the origin is substantially increased. 
Fitting instead its value to the mass difference between 
J/ip and rjc, we find Ccc/rri?c — ^^'^ MeV. Assuming that 
the corresponding attraction between charmed diquark is 
three times smaller than that between the charm quark- 
antiquark pair as in the case of light quarks, we have 
Ccc/rn^ — 39 MeV. We could introduce additional mass 
dependence in Cs and in Cm by fitting the mass difFer- 
ences in the strange, charm and bottom hadrons from 
Table I and II, respectively. However, these introduce 
only minor changes in the analysis to foUow, and there- 
fore we will just use the mass independent C/f 's obtained 
above. 



B. charmed tetraquark mesons 

Using above parameters, we argue in this subsection 
why the doubly charmed tetraquark meson might be sta- 
ble. Let us consider a tetraquark meson Tq-^q^ that is 
made up of udqiq2- The reason we start with the ud 
diquark is that for a diquark the strongest attraction is 
expected when the two quarks are light, and their total 
color, flavor and spin are all in the antisymmetric states. 
Therefore, if there is any stable configuration, it must 
involve a scalar ud diquark. We then add two antiquarks 
in the relative s-wave statě and look for a stable config- 
uration. 



The stability of Tq-^q^ depends on whether it is ener- 
getically favorable against recombining into two mesons 
of uqi and dq2- As we have discussed previously, the at- 
traction Cm between a quark-antiquark pair is stronger 
than Cb in a diquark. This means that when both qi 
and q2 are light, the two- meson states would be energet- 
ically much more favorable, and Tq-^q^ will not be stable. 
However, when qi and q2 become heavy, the attraction in 
the quark-antiquark pair in the meson decreases, while in 
^9i<32 the attraction in the ud diquark remains the same 
and the interaction in the qiq2 decreases substantially. 
Therefore, the tetraquark statě could become stable. A 
simplification in working with a spin zero ud diquark in 
Tq-^q^ is that there is no spin-spin interaction between the 
ud diquark and qi or q2, and it is sufhcient to only esti- 
mate the attractions inside the diquark or anti-diquark. 
If qi and q2 are identical quarks, then their total spin 
has to be zero, because their color combination is an- 
tisymmetric in the present configuration. This means 
that their total spin has to be 1, which is a repulsive 
combination. However, the repulsion becomes smaller 
when quark masses become heavy. Moreover, the quan- 
tum number of Tq-^q^ has to be 1"*", so that it can not 
decay into two pseudoscalar mesons. The threshold for 
its decay is then the masses of the vector and pseudo 
scalar mesons. 



udqiq2 (spinal) 

3 Cb i 1 Cb 

4 mj 4 m^^ 


uqi (spinal) 

1 Cm 


dq2 (spin=0) 

3 Cm 


udqiq2 

—uqi — uq2 


4 m^tmg^ 


4 m-u^mq-^ 


udšs 


K* 


K 




-127 


92 


-285 


63 


udcc 


D* 


D 




-143 


31 


-95 


-79 


udhh 


B* 


B 




-145 


10 


-30 


-124 



TABLE IIL Tetraquark mesons udqiq2 with spin equal tol 
for qi — q2, where qi, q2 = s, c and b. Units are in MeV. 

Table IIIII shows the mass difference between a 
tetraquark meson with identical diquarks and the sum 
of vector and pseudo scalar meson masses due to the 
color-spin interaction of Eq. ([l]) with the Ch parameters 
determined previously. As expected, the mass difference 
decreases as qi and q2 become heavy, and the tetraquark 
mesons Tec and Tbb with c or b quarks are bound. Al- 
though our result is based on a very crude estimate, es- 
sentially the same result has been obtained in the fuU 
constituent quark model calculation [S] |TF; and the QCD 
sum-rule calculation [16 . 

For qi and q2 of different fiavors, their total spin could 
be either zero or one. The quantum number of the 
tetraquark meson could then be either 0+ or 1+. Ta- 
ble IIVI and Table |V] show the mass difFerences in such 
cases. As in the previous case, bound tetraquark mesons 
with čb could exist. 
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udq\q2 (spin=0) 


ngi (spin=0) 


ag2 (spin=0) 


udqiq2 


ududq 


uud 


dq2 


ududq — uud — dq 


3 Cb 3 Cb 


3 Cm 


3 Om 


~uqi — uq2 


2(-!£t) +^^^=1 
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3 C M 
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4 m^tmg2 








udsč 




D 




ududs 


N 






-162 


-285 


-95 


218 




-145 


-286 


141-KASl=i 


udsč 


K 






ududč 


N 


D 




-150 


-285 


-30 


165 




-145 


-95 


-50+Aí;í,=i 


udcb 
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B 




ududb 


N 


B 




-146 


-95 


-30 


-21 




-145 


-30 


-114+A£l=i 



TABLE IV: Tetraquark mesons udq\q2 with spin equal to O 
for q\^ q2- qi, q2 = s, c and b. Units are in MeV. 



udq\q2 (spinal) 

3 Cb 1 1 Cb 

4 4 rriq-^ 


uq\ (spinal) 
1 Cm 

4 irLumq-^ 


dg2 (spin=0) 

3 Cm 

4 m^rriq^ 


udqiq2 
~uqi — uq2 


itdsč 
-139 


K* 
95 


D 

-95 


-139 


_D* 
31 


A' 
-285 


114 


udsb 
-143 


K* 
95 


B 
-30 


-208 


B* 
10 


K 

-285 


132 


udčb 
-144 


_D* 
31 


B 
-30 


-145 


B* 
10 


D 

-95 


-59 



TABLE V: Tetraquark mesons udqiq2 with spin equal to 1 
for qi 7^ g2, where qi, q2 = s, c and b. Units are in MeV. 

C. charmed pentaquark baryons 

Similar observations can be made for heavy pentaquark 
baryons. Many constituent quark model calculations 
show that the 0+ [T7], if it exists at all, can not be 
explained as a bound statě of ududš constituent quarks 
|18j . This is due to the strong attraction between the 
š and the hght quark, so that it is energetically much 
more favorable for ududš to form a meson and a baryon. 
The attraction to form a meson becomes smaller if the 
š is replaced by either a č or 6. Full constituent quark 
model calculations [121 HH] indeed find a possible stable 
heavy pentaquark baryon. A likely pentaquark struc- 
ture would be that suggested in Ref. [^I] with the two 
scalar diquark ud combined into an L = 1 and color an- 
tisymmetric statě. The excitation energy of a diquark 
in a L = 1 statě, AE^^i, can be estimated by approxi- 
mating the charmed baryon as a sum of a charm quark 
and a diquark, because the interaction between thcm is 
small in the heavy quark limit. Attributing the mass dif- 
ference between the parity doublet partners of the posi- 



TABLE VL Charm and bottom pentaquark baryons 
Qq{ududq) {q = c and b) with spin equal to 1. AE^^^ — 309 
MeV is en excitation energy of two diquarks with relative an- 
gular momentům L = 1. Units are in MeV. 



tive parity A+ (2286 MeV) and the negative parity A*+ 
(2595 MeV) to the L = 1 excitation of the diquark, as 
the heavy charm quark would act as the center of mass, 
leads to AEl=i = 309 MeV. Applying this L = 1 exci- 
tation energy to the relative excitation of two diquarks, 
we find that while a strange pentaquark baryon is very 
unlikely to exist, heavy pentaquark baryons &c and Of, 



could be closer to the threshold as shown in Table VI 
consistent with the full constituent quark model calcu- 
lation [m [ini US and the QCD sum-rules study [13] 
in which a possible stable heavy pentaquark baryon has 
been found. 
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-145 
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TABLE VII: Charm- and bottom-strange pentaquark baryons 
Qqsiudusq) {q — c and b) with spin equal to 0. Units are in 
MeV. 

For a pair of ud and us diquarks in Qcs{udusc), they 
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do not have to be in the L = 1 statě, and hence there is 
no additional contribution from the orbital energy |24j . 
The result from our simple estimates are given in table 
VII Previous experiments [251 12S] have tried to search 



In this model, the Tec number is given by 



for this pentaquark baryon assuming that it is bound and 
has a lifetime similar to that of Ds ■ The experiment could 
only determine an upper bound of greater than 0.02 for 
its production cross section relative to that for the Ds, 
which is larger than typical theoretical estimates. From 
simple application of statistical hadronization model, the 
number of &čs to that of Dg is roughly exp(— (mog^ — 
m£)^)/T) ~ exp(— 5) = 0.007, assuming a hadronization 
temperature of T = 200 MeV. This is smaller than the 
experimental upper bound and therefore further search 
is essential. 



III. PRODUCTION OF CHARMED EXOTICS 
IN RELATIVISTIC HEAVY ION COLLISIONS 



As discussed above, tetraquark mesons and pentaquark 
baryons are more likely to exist in the heavy quark sec- 
tor such as the Tec, Tet, Thb, Qcs, and Qc- It is, however, 
very unlikely that they can be observed in B decays or 
elementary processes, as the favorable exotics involve two 
heavy quarks. However, the abundance of heavy quarks 
is significantly enhanced in ultrarelativistic heavy ion col- 
lisions, e.g. in a centrál coUision at the LHC, more than 
20 cč pairs are expected to be produced in one unit of 
midrapidity. Therefore, while a heavy quark produced in 
an elementary process will most likely find a heavy anti- 
quark instead of a heavy quark, the probability to find a 
heavy antiquark or a heavy quark in a heavy ion coUision 
is similar. The probability to form a Tec compared to a 
J/ip thus will only be suppressed by the additional sta- 
tistical factor coming from combining an additional ud 
diquark. 



A. charmed tetraquark mesons 

The number of heavy tetraquark mesons produced 
from the quark-gluon plasma formed in relativistic heavy 
ion coUisions can be estimated in the coalescence model 
[27j , which has been shown to describe very well the pion 
and proton transverse momentům spectra at intermedi- 
ate momenta [2SJ [53] as well as at low momenta if res- 
onances are included [30], and the yield and transverse 
momentům spectra of phi meson and Omega baryon [31 
as well as the charmed meson [32]. We employ the for- 
mula that was previously used to calculate the yields 
of tetraquark Dsj{2317) meson [33 and pentaquark 0+ 
baryon [T2] at RHIC to study Tec production in centrál 
Au+Au coUisions at RHIC and Pb+Pb coUisions at LHC. 



coal 



x/T„(a^i--a^4;Pi--P4)- 



(2) 



In the above, the color-spin-isospin factor g^^^ = 3 x 
1/3^ X 1/2* = 1/432 is the color-spin-isospin factor for 
the four quarks to form a hadron of the quantum num- 
ber of the tetraquark meson and da denotes an element 
of a space-like hypersurface at hadronization. Assuming 
Bjorken correlation y = rj between the space-time rapid- 
ity r] and the momentum-energy rapidity y and neglect- 
ing the transverse flow as well as using the non-relativistic 
approximation, we obtain the foUowing expression for the 
number of Tec produced from quark coalescence: 



1 N-c N-cN^Ng A i^naf)^/^ 

2 l\Vcil + 2^,,Tc<Jfy 



432 



where Tc = 170 MeV is the critical temperature and Vc 
is the fireball volume at hadronization, which is about 
1,000 fm^ in centrál Au-|-Au coUisions at s}/^ = 200 GeV 
[53] and about 2,700 fm'' in centrál Pb+Pb coUisions at 

1/2 

^NN ~ '^■'^ TeV [11. . The quark numbers at hadroniza- 
tion are denoted by and Nd for light quarks and Nc 
and Nč for heavy quarks. Their values are taken to be 
N^ = Nd = 245 [37 and 662 TT] as well as Nc = Ns = 5 
and 20 in centrál RHIC and LHC coUisions, respectively, 
all in one unit of midrapidity. The charm quark numbers 
are based on initial hard scattering of nucleons in the col- 
liding nuclei [II1[33|- In obtaining Eq. we have used 
the quark momentům distribution function 

f,{x,p) = 6<5(77 - y)exp(-(m2 +p2,)i/2/Te) (4) 

and the tetraquark meson Wigner distribution function 



fTj^;p) 



8^ exp 



(5) 



1=1 



where the relative coordinates and momenta are 
related to the quark coordinates and momenta by 
the Jacobian transformations defined in Eqs. (7) and (8) 
of Ref. [35]. The width parameter CTí in the Wigner 
function is related to the oscillator frequency lo hy ai — 
í/y/JiiUj with the reduced masses /íí defined in Eq. (9) of 
Ref. [33]. 

In Fig. [T] we show the numbers of Tec produced at 
RHIC and LHC as functions of the oscillator frequency. 
Because of the larger abundance of charm quarks at LHC 
than at RHIC, the number of Tec produced at LHC is 
more than an order of magnitude larger than that pro- 
duced at RHIC. For the oscillator frequency lu — 0.3 
GeV, determined from the size (r|) ,)ch « 0.124 fm^ of 
the D+{cš) meson based on the light-front quark model 
[33] , the number of Tec produced at RHIC and LHC is 
about 5.5 X 10~^ and 9.0 x 10~^, respectively. 
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FIG. 1: Numbers of T^c produced at RHIC and LHC as func- 
tions of the oscillator frequency used for the quark wave func- 
tions in Tec 



It is of interest to compare the predicted number of 
Tec mesons from the coalescence model with that from 
the statistical model. In this model the number of T^c 
mesons produced at hadronization is given by |33j : 



6cs produced at RHIC and LHC are about 1.2 x 10~^ 
and 7.9 x 10"'', respectively. 

Since the predicted numbers of Ds meson from the 
coalescence model are about 5.3 x lO^'^ at RHIC and 0.58 
at LHC, the estimated ratio of numbers of Qcs and is 
about 2.3 X 10"^ at RHIC and LHC. This is consistent 
with the Boltzmann factor due to the uud component 
in &CS- In fact, extracting sč component from uudsč, 
the remaining uud N component has a Boltzmann 
factor e"™"/^ ~ 4.0 x 10"^ with T = 170 MeV. Similar 
estimate also works for the case of the A and D in Qcs- 
Using the value 0.16 and 1.1 for the D meson numbers 
at RHIC and LHC, respectively, the calculated ratio of 
numbers between Qcs and D is about 0.74 x lO"'^, while 
the uds component has a Boltzmann factor of er™-'^!^ ~ 

1.4 X 10-3. 

In the statistical model, the yield of Qcs is given by a 
formula similar to Eq. (|6| except the power in the charm 
fugacity parameter 7c. Since there is only one charm 
quark in Qcs-, the yield is only proportional to ^c- Using 
same parameters for evaluating the yield of Tec, we obtain 

4.5 X 10"^ and 2.7 x 10"^ for Qcs produced in centrál 
Au+Au coUisions at RHIC and centrál Pb+Pb coUisions 
at LHC, respectively. These values are again significantly 
larger than those predicted from the coalescence model. 



^stat 



(27r)2 



dmTmj,e /g I — — ) (6) 



where Vh and (3h are the volume and radiál flow velocity 
of formed hadronic matter, and 7^7 is the fugacity param- 
eter for ensuring that the number of charmed hadrons 
produced statistically at hadronization is same as the 
number of charm quarks in the quark-gluon plasma. 
With Vh ~ 1,908 fm^ Ty = 175 MeV, = 0.3c, 
and the charm fugacity 7c « 8.4 [33], we obtain Nt^^ ~ 
7.5 X 10"'' in centrál Au+Au coUisions at RHIC. The 
yield of Tec increases to 8.6 x 10"^ in centrál Pb+Pb 
coUisions at LHC where we have used Vh ~ 5,220 fm'^, 
Th = 175 MeV, Ph = 0.47c, and the charm fugacity 
7c ~ 16.3 [11]. Compared to those from the coalescence 
model, predictions from the statistical model are almost 
two orders of magnitude larger. 



B. charmed pentaquark baryons 

For the yield of pentaquark baryon Qcs{udusc), the 
coalescence model gives 



1 NsN^NuNd -rV {^^cjW^ 



Ayc(l + 2/i.Tca2)' 



Using again the oscillator frequency oj — Q.?> GeV and 
taking the anti-strange quark numbers to be 150 |33j and 
405 [TP at RHIC and LHC, respectively, the numbers of 



IV. DECAY MODES OF CHARMED EXOTICS 

In this section, we discuss the observable decay modes 
of the tetraquark Tec and the pentaquark Qcs- As we 
have discussed already. Tec is most Ukely a stable statě, 
since its mass is below the threshold of D*D. To be 
more generál, we consider nevertheless both cases where 
the mass of Tec is above or below the threshold, and dis- 
cuss in each case possible decay modes that can be real- 
istically detected in experiments with good performance. 
For the Tec above the threshold of D*D, it can decay 
to D*^ via a strong process [40 . For the Tec below 
the threshold of D*D and above DDtt, the decay chan- 
nel to D*^ is energetically forbidden, but the D*^ 
component in Tec can decay through a strong process, 
leading to the finál decay mode Ď^Ď^'k~ . On the other 
hand, when Tec is below the threshold of DDtt, the de- 
cay channel of D*~ is closed and only the weak decay 
of Z)° component in Tec is allowed via K~^Tr~ or 

/í+7r+7r"7r". Therefore, Tec would be detected by the 
decay modes D*~K'^tt~ and D*~ K^Tr^^Tr'^' . The last 
two decay patterns would most likely happen since the 
binding energy of Tec is estimated to be about 80 MeV 
as shown previously, which is sufficiently larger than the 
mass difference (about 6 MeV) between D*^ and Ď^ii^ . 
Below the threshold of DDtt, it may be also interesting 
to see the decay of D*~ component in Tec- Consider- 
ing that D*~ component contains a quantum number of 
Ď'^TT~, and decays into K~^tt~ and iř+7r+7r"7r", we 



may observe the D^K^ 



and D°K 



TT TT TT TT 



de- 
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cays. 

Among the weak decays below the threshold of DDtt, 
the decay of the l)° componcnt in Tec can be dis- 
tinguished from that of the D*~ component. The 
formcr has the correlations {K~^TT~){K^Tr~)TT~ and 
(íí+7r+7r+7r^)(iř+7r^)7r^, and the latter has the correla- 
tions {K^Tr^){K^TT^Tr~) and (if +7r~)(íí+7r+7r~7r~7r~), 
where brackets denote correlated particles. However, the 
D^D^TT~ statě, which would appear in Tec in the lat- 
ter process, contains six quarks, hence further analysis is 
needed to discuss its stability. 

The pentaquark @cs also has interesting decay pat- 
terns. As can be seen in table IVIIl the mass of Qr^ could 



be slightly above the A£>° threshold, in which case its 
lifetime will be shorter than that of Dg- Then, the only 
possible way to look for it is from the hadronic decay 
to A + hnal states. Since ALICE will be able to re- 
construct the Ď'^ through its hadronic decay, it will be 
an excellent opportunity to search for Qcs- Considering 
more generál cases, and assuming Qcs to be above the 
threshold of NDg, the Qcs can decay into pD~ and Aí)° 
or AD~ via the strong process. Although the SZ? chan- 



nel is also a possible decay mode, it is more difíicult to 
detect as compared to NDg and AD. When the mass 
of Qcs is below the ND^ and above the AD threshold, 
it decays only to AĎ^ or AD^ . On the other hand, be- 
low the threshold of AD, the hadronic decay channels 
are closed and only weak decays are possible. In this 
case, the lifetime of Qcs will depend on the lifetime of 
the different components inside the Qcsi such as the A, 
and D~ , whose lifetimes are respectively 2.6 x 10~^°, 
0.41 X 10~^^ and 1.0 x 10~^^ seconds. Therefore, once 
the Qcs is formed as a deeply bound statě, it will decay 
by the weak process of Ď° or D~ . Consequently, pos- 
sible finál states would be AK^tt^, Aif +7r+7r~7r~ and 
AK+TT^n^. 

Since the lifetimes of Tec and Qcs are in the order of 
10^^^ seconds, their decays occur outside the coUision re- 
gion and they are thus identifiable by vertex reconstruc- 
tion. Therefore, Tec and Qcs would be identified clearly 
in experiments if they exist. We summarize our results 
on possible decay modes of Tec and Qcs in Tables VIII 
andHXl 



TABLE VIII: Possible decay modes of T^c- In the bottom col- 
umn, we would observe correlations {K'^TT~){K'^n~)n~ and 
{K^-K^Tr'^TT~){K'^-K~)TT~ In the finál states. See the text for 
details. 



threshold 


decay mode 


lifetime 


AÍT,, > Md' + Md 


D*-D° 


hadronic decay 


2Md + < Mt,, < Md* + Md 




hadronic decay 


Mt,, < 2Md + 


D*-K+Tv', D*'K+TV+TT-TV- 


0.41 X 10"^^ sec. 



TABLE IX: Possible decay modes of 6c 



threshold 



decay mode 



lifetime 



Me,, > Mn + Md, 
Ma + Md < Me„ < Mn + Md, 

Me„ < Ma + Md 



pD7 
AD- 

AK + TV-, AK+TV + TV-TT- 
ATf+TT^TT" 



hadronic decay 
hadronic decay 
hadronic decay 
0.41 X 10"^^ sec. 
1.0 X 10"^^ sec. 



Lastly, we comment on the possibility to measure dou- 
bly charmed baryons in heavy ion coUisions. The doubly 
charmed baryon 2++ have been observed by the SELEX 
CoUaboration in the A+if ~7r+ and in the pD~^K~ decay 
modes with a mass of (3518.7 ± 1.7) MeV [35l [3^. The 
same coUaboration has also successfuUy measured S+ in 



the A+íí 7r"'"7r+ decay mode with a mass of 3460 MeV 
[37] . On the other hand, attempts by the FOCUS Col- 
laboration in the photoproduction experiment, and by 
the BAB AR in e+e" annihilation experiments [311 ISS] 
have SO far failed to establish the existence of the dou- 
bly charmed baryons. Hence, it is an interesting problém 
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to search for the doubly charmed baryons in heavy ion 
coUisions. Using the coalescence model, we find that the 
number of S+ produced are 1.9 x 10^^ at RHIC, and 
3.2 X 10~^ at LHC. Therefore, we will be able to reahsti- 
cally measure S+ and through their decay vertices 
to A+iír^7r+ and pD^K^ , and to A+/í^7r+7r+, respec- 
tively. 

V. SUMMARY 

Based on the consideration of the color-spin interaction 
between diquarks, which describes reasonably the mass 
splittings between many hadrons and their spin fiipped 
partners, we have shown that tetraquark mesons and 
pentaquark baryons that consist of two charmed quarks 
could be bound. Using the quark coalescence model, 
their yields in heavy ion coUisions at both RHIC and 
LHC are estimated. Because of the expected large charm 
quark number in centrál Pb+Pb coUisions at LHC, the 
abundances of the tetraquark meson Tec and pentaquark 
baryon Qcs are about 10^^ and 10^'^, respectively. We 
have also discussed their decay modes to illustrate how 
they can be identified in heavy ion coUisions. In our stud- 
ies, we have not taken into account the hadronic effect 
on the abundance of these charmed exotics, as hadronic 
reactions that affect their annihilation and production 



are unknown. Since the yields of T^c and Qcs from the 
coalescence model is significantly smaller than those ex- 
pected from the statistical hadronization model, includ- 
ing the hadronic effect is expected to increase their yields 
substantially and reduces the differences from the predic- 
tions from the quark coalescence model and the statistical 
hadronization model. Also, charmed hadrons would be 
more abundantly produced, particularly the Tec, if charm 
quarks are produced from the QGP formed in these col- 
lisions. We also comment on the possible measurement 
of doubly charmed baryons in heavy ion coUisions, and 
the estimated numbers are 1.9 x 10~^ and 3.2 x 10"'' at 
RHIC and LHC, respectively. We thus expect that the 
open and hidden charmed hadron physics will be an in- 
teresting subject in the forthcoming heavy ion coUision 
experiments. 
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